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ABSTRACT 

If supermassive black holes (SMBHs) are the energy sources that power quasars and active galactic 
nuclei, then QSO SDSS 1148+5251, the quasar with the highest redshift (zqso = 6.43), hosts a 
supermassive black hole formed within ~ 0.9 Gyr after the Big Bang. This requirement places 
constraints on the cosmological formation of SMBHs, believed to grow from smaller initial seeds by a 
combination of accretion and mergers. We focus on gas accretion onto seeds produced by the collapse 
of Pop III stars at high redshift. We incorporate the results of recent relativistic, MHD simulations 
of disk accretion onto Kerr black holes to track the coupled evolution of the masses and spins of the 
holes. We allow for an additional amplification of ~ 10 4 in the mass of a typical seed due to mergers, 
consistent with recent Monte Carlo simulations of hierarchical mergers of cold, dark matter halos 
containing black hole seeds. We find that the growth of Pop III black hole remnants to ~ 10 9 M© 
by zqso ^ 6.43 favors MHD accretion disks over standard thin disks. MHD disks tend to drive the 
holes to a submaximal equilibrium spin rate a/M ~ 0.95 and radiation efficiency 6m ~ 0.2, while 
standard thin disks drive them to maximal spin (a/M = 1) and efficiency (cm = 0.42). This small 
difference in efficiency results in a huge difference in mass amplification by accretion at the Eddington 
limit. The MHD equilibrium efficiency is consistent with the observed ratio of the QSO plus AGN 
luminosity density to the local SMBH mass density. Our prototype analysis is designed to stimulate 
the incorporation of results from relativistic stellar collapse and MHD accretion simulations in future 
Monte Carlo simulations of hierarchical structure formation to better determine the cosmological role 
of SMBHs and their mass and spin distributions. 

Subject headings: black hole physics — accretion — MHD — cosmology — quasars 



1. INTRODUCTION 

There is substantial evidence that supermassive black 
holes (SMBHs) with masses in the range 10 6 - 10 10 M Q 
exist and are the engines that power active galactic nu- 
clei (AGNs) and quasars (Rees 1984, 1998, 2001; Mac- 
chetto 1999). There is also ample evidence that SMBHs 
reside at the centers of many, and perhaps most, galax- 
ies (Richstone et al. 1998; Ho 1999), including the Milky 
Way (Genzel et al. 1997; Ghez et al. 2000, 2003; Schodel 
et al. 2002). 

The highest redshift of a quasar discovered to date is 
2QSO = 6.43, corresponding to QSO SDSS 1148+5251 
(Fan et al. 2003). Accordingly, if they are the en- 
ergy sources in quasars (QSOs), the first supermassive 
black holes must have formed prior to zqso = 6.43, or 
within t — 0.87 Gyr after the Big Bang in the concor- 
dance ACDM cosmological model. This requirement sets 
a significant constraint on black hole seed formation and 
growth mechanisms in the early universe. Once formed, 
black holes grow by a combination of mergers and gas 
accretion. An important clue to the growth process is 
provided by the ratio R of the observed QSO plus AGN 
luminosity density to the local SMBH mass density, since 
R is related to the mean radiative efficiency tu of accre- 
tion onto black holes (Soltan 1982): cm > R- Recent 
measurements suggest 0.1 < R < 0.2 (Yu & Tremaine 
2002; Elvis, Risaliti & Zamorani 2002), a range consis- 
tent with disk accretion onto black holes. Another clue 
is provided the estimated ratios cl of the bolometric- 



to-Eddington luminosities of the broad-line quasars in 
a Sloan Digital Sky Survey sample of 12,698 quasars in 
the redshift interval 0.1 < z < 2.1. This survey sup- 
ports the value « 1 as a physical upper limit (McLure 
& Dunlop 2004). Barring an extraordinary coincidence, 
the range of values inferred from observations for both 
cm and e l suggest that a significant fraction of the mass 
of SMBHs is acquired by gas accretion. Together, these 
two parameters control the rate of growth of black holes 
by accretion and are crucial in determining whether or 
not initial seed black holes formed at high redshift have 
sufficient time to grow to SMBHs to explain quasars and 
AGNs. 

The more massive the initial seed, the less time is re- 
quired for it to grow to SMBH scale and the easier it is to 
have a SMBH in place by z > 6.43. One possible progeni- 
tor that readily produces a SMBH is a supermassive star 
(SMS) with M > 10 3 M Q (see, e.g., Shapiro 2004 for 
a recent review and references). SMSs can form when 
gaseous structures build up sufficient radiation pressure 
to inhibit fragmentation and prevent normal star forma- 
tion; plausible cosmological scenarios have been proposed 
that can lead to this occurence (Gnedin 2001; Bromm 
& Loeb 2003). SMSs supported by radiation pressure 
will evolve in a quasi-stationary manner to the point of 
onset of dynamical collapse due to a general relativis- 
tic radial instability (Chandrasekhar 1964a, b; Feynman, 
unpublished, as quoted in Fowler 1964). The collapse 
of a nonrotating, marginally unstable, spherical SMS of 
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mass M leads directly to the formation of a nonrotating 
Schwarzschild black hole of the same total mass (Shapiro 
& Tcukolsky 1979). But like most stars formed in na- 
ture, SMSs will be rotating. In fact, in the event that 
viscosity and/or turbulent magnetic fields are present to 
drive these stars to uniform rotation, they are likely to 
be maximally rotating (i.e. at the mass-shedding limit) 
by the time they reach the onset of collapse (Baumgarte 
& Shapiro 1999). Recent relativistic hydrodynamic sim- 
ulations have shown that unstable, maximally rotating 
SMSs of arbitrary mass M inevitably collapse to SMBHs 
of mass ~ 0.9M and spin parameter a/M ~ 0.75; the 
rest of the mass goes into an ambient disk about the 
hole (Shibata & Shapiro 2002: Shapiro & Shibata 2002; 
Shapiro 2004). 

But the fact remains that SMSs have yet to be ob- 
served and there is no concrete evidence that they actu- 
ally form in the early universe. Moreover, simulations of 
cosmological structure formation performed to date indi- 
cate that the first generation of stars are more likely to be 
zero-metallicity Pop III stars in the range 10 2 — 10 3 M Q 
(Bromm, Coppi & Larson 1999, 2002; Abel, Bryan, & 
Norman 2000, 2002; but note that Norman 2004 re- 
ports preliminary indications that supermassive stellar 
objects with M > 10 4 M Q may be forming as second 
generation stars at 10 < z < 15 in his latest simula- 
tions). So the most conservative hypothesis is that the 
seed black holes that later grow to become SMBHs orig- 
inate from the collapse of Pop III stars (Madau & Rees 
2001) and not SMSs, and this is the ansatz we shall ex- 
plore here. Newtonian simulations suggest that Pop III 
stars with masses in the range M ~ 60 — 140M Q and 
M > 260M Q collapse directly to black holes, while stars 
with M ~ 140 — 260M© undergo explosive annihilation 
via pair-creation processes (Heger et al. 2003). The up- 
per limit to the mass of a Pop III star is set by accretion 
over a stellar lifetime, yielding M < 600M Q (Onukai & 
Palla 2003; Yoshida et al. 2003). Abel et al. (2002) 
argue that first generation stars significantly larger than 
100M Q are likely to explode before they have time to 
accrete to larger masses. The Pistol Star is an example 
of an existing star believed to have a mass > 200M Q , 
although it has high metallicity (Figer et al. 1998). Very 
massive stars are dominated by thermal radiation pres- 
sure, so the catastrophic collapse of those that do not 
explode will be hydrodynamically similar to the collapse 
of SMSs, producing black holes with masses comparable 
to those of their progenitors. The higher the mass of 
the black hole seed, and the earlier it forms in the uni- 
verse, the easier it is for it to grow to a SMBH, hence 
the more conservative will be any contraints imposed on 
the cosmological black hole growth rate by the existence 
of a SMBH by zqso = 6.43. Accordingly, we take the 
highest range of plausible values for the masses of black 
hole seeds to establish the most conservative (robust) 
constraints, adopting 100 < M/M & < 600 for the range 
of black hole seeds arising from the collapse of Pop III 
stars at z < 40. 

The most detailed studies of SMBH formation to date 
involve detailed, Monte Carlo simulations that follow 
the cosmological growth of a distribution of black hole 
seeds by a combination of discrete, stochastic mergers 
as well as gas accretion (see, e.g., Haehnelt & Kauff- 
mann 2000; Volonteri, Haardt & Madau 2003; Brom- 



ley et al. 2004; Haehnelt 2004; Haiman 2004; Yoo & 
Miraldi-Escudc 2004; and references therein). Both pro- 
cesses are assumed to take place in the context of the 
cold dark matter (CDM) model, where dark matter halos 
merge hierarchically, and black holes are assumed to set- 
tle, merge and accrete in their gaseous centers. Typically, 
the stellar dynamical processes that lead to mergers, as 
well as the hydrodynamical processes that fuel accretion, 
are modeled in these analyses by implementing simple, 
physically plausible, rules rather than by detailed inte- 
grations of the governing dynamical and hydrodynamical 
equations of motion. Performing such "first principles" 
integrations would prove prohibitive in this context. 

In this paper we focus on SMBH cosmological growth 
by accretion. We identify and explore the main variables 
that govern this process and ultimately influence the 
outcome of detailed Monte Carlo simulations that track 
SMBH growth. We incorporate some of the most recent 
findings of relativistic magnetohydrodynamical (MHD) 
simulations of gas accretion onto black holes and ex- 
plore their cosmological implications for SMBH evolu- 
tion. We show how the evolution and amplification of 
black hole mass by accretion is intimately tied to the 
evolution of black hole spin, probing some of the impli- 
cations of our earlier survey and analysis of black hole 
spin evolution (Gammie, Shapiro & McKinney 2004) for 
cosmology. We demonstrate how in principle the very 
existence of a quasar at redshift zqso = 6.43 can help 
constrain the formation epoch and/or size of black hole 
seeds and select among competing models of accretion. 

Specifically, recent relativistic MHD simulations pre- 
dict the radiation efficiency cm as a function of the BH 
spin parameter a/M; they also predict the spin-up rate 
as a function of a/M. Monte Carlo simulations that de- 
termine cosmological SMBH growth must integrate cou- 
pled mass and spin evolution prescriptions versus time 
for each accreting BH to produce reliable growth his- 
tories and final SMBH masses and spins. This paper 
provides prototype integrations of the coupled evolution 
equations, focussing on the history of the progenitor of 
SDSS 1148+5251 and demonstrating the significant dif- 
ferences in the outcomes for standard thin disk accretion 
models versus recent MHD models. 

Our discussion is simplified and illustrative at best: 
much of the input physics involving black hole seed for- 
mation, accretion flows, and mergers is still being devel- 
oped. Our main goal is to isolate some of the underlying 
local physical issues and parameters pertaining to accre- 
tion to better understand their role in determining the 
global outcome of the cosmological Monte Carlo simu- 
lations of SMBH build-up during hierarchical structure 
formation. Excellent overviews of the input physics have 
appeared elsewhere (see, e.g. Haiman & Quatacrt 2004, 
and references to earlier work). Also, earlier treatments 
have considered some of the constraints on cosmologi- 
cal growth imposed by the recent discoveries of lumi- 
nous quasars at high redshift (e.g., Haiman & Loeb 2001; 
Haiman 2004; Yoo & Miraldi-Escude 2004). But here we 
specifically want to illustrate in the simplest fashion how 
the most recent findings related to relativistic, MHD ac- 
cretion flows onto spinning black holes have important 
implications for evolutionary models of the growth of 
SMBHs in the early universe. We emphasize by con- 
crete example the point made in Gammie et al. (2003) 
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that tracking the spin as well as the mass of a black holes 
is necessary to determine its growth (see also Volonteri 
et al. 2004). We also show that the value of accretion 
radiation efficiency, ~ O.f , adopted in many Monte Carlo 
simulations may not be entirely consistent with the latest 
MHD accretion disk modeling. Determining this param- 
eter, on which the outcome of cosmological simulations of 
SMBH growth depends very sensitively (exponentially!), 
is coupled to the spin evolution of the hole; both may 
now be within our grasp via detailed relativistic MHD 
simulations of BH accretion. 

The calculations performed here are prototypical only; 
our main aim is to motivate the incorporation of these 
parameters in more detailed Monte Carlo studies and 
thereby sharpen some of the rules that enter these sim- 
ulations. We also hope to provide those members of 
the computational relativistic MHD community who may 
not be SMBH Monte Carlo cognoscenti a simple means 
of extracting the essence of the Monte Carlo simulations, 
particularly the evolutionary consequences of different 
accretion models. 

In Section 2 we set out the basic equations that de- 
scribe the coupled evolution of black hole mass and spin 
by accretion. We also summarize in this section the 
relations we require from the concordance ACDM cos- 
mological model. In Section 3 we integrate the coupled 
evolution equations to track the evolution of black hole 
mass and spin as functions of time. In Section 4 we ap- 
ply the results to the cosmological problem and identify 
some constraints imposed by the existence of a SMBH at 
z > 6.43. In Section 5 we summarize briefly and discuss 
some caveats and areas for futher study. 

2. BASIC EQUATIONS 

Here we assemble the fundamental black hole accretion 
evolution equations and review the underlying assump- 
tions on which they are based. We then specify a back- 
ground cosmological model in the which the growth of 
the black hole to supermassive size occurs. 

2.1. Black Hole Mass and Spin Evolution 

Define em, the efficiency of conversion of rest-mass en- 
ergy to luminous energy by accretion onto a black hole 
of mass M, according to 



e M = L/M c 2 



(1) 



where Mq is the rate of rest-mass accretion and L is the 
luminosity. Define el, the efficiency of accretion lumi- 
nosity, according to 



e L = L/L E , 



(2) 



where Le is the Eddington luminosity, given by 

L E = ^ « 1.3 x 10 4 VM 8 erg s -i. (3 ) 

Here we assume that the accretion is dominated by nor- 
mal, baryonic matter and ignore any contribution of col- 
lisionless or self- interacting dark matter (but see, e.g., 
Ostriker 2000 and Balberg & Shapiro 2002 for alternative 
scenarios). In particular, we assume that the accreting 
gas consists of fully ionized atoms and that the princi- 
pal opacity source is Thomson scattering. The quantity 
H e is the mean molecular weight per electron and m p is 



the proton mass. The black hole growth rate must ac- 
count for the loss of accretion mass-energy in the form 
of outgoing radiation according to 



W = (1 - e M )M . 



(4) 



Combining eqns. Q-Q we obtain the black hole growth 
rate, 

dM _ e L (l - cm) M 
dt em t ' 

where r is the characteristic accretion timescale, 
Mc 2 



(5) 



0.45/j- 1 Gyr, 



(6) 



(8) 



and is independent of M. 

The mass accretion efficiency em is typically a func- 
tion of the black hole spin parameter a/M = J/M 2 . It 
changes with time as the spin evolves. It is convenient to 
express the spin evolution in terms of the nondimensional 
quantity s = s(a/M), defined by 

g _ d(a/M) M 

dt M 

Inserting eqns. J3J and © into eqn. {JJ yields the evolu- 
tion equation for the black hole spin, 

d{a/M) _ e L s 

dt cm t 

In general, eqns. JSJ and JH} must be integrated simulta- 
neously to determine the mass and spin evolution of the 
black hole. 

Determining ej\j(a/M) and s(a/M), which are needed 
to integrate eqns. JSJ and JSJ), requires a gas dynam- 
ical model for black hole accretion. We shall assume 
that the gas has sufficient angular momentum to form a 
disk about the hole and consider two different accretion 
disk models: (1) a standard, relativistic, Keplerian "thin 
disk" with "no-torque boundary conditions" at the inner- 
most stable circular orbit (ISCO) (Pringle & Rees 1972; 
Novikov & Thorne 1973; see Shapiro & Teukolsky 1983 
for review and references) and (2) a relativistic, MHD ac- 
cretion disk that accounts for the presence of a frozen-in 
magnetic field in a perfectly conducting plasma (De Vil- 
liers & Hawley 2003; De Villiers, Hawley & Krolick 2003; 
Gammic, Shapiro & McKinney 2003; McKinney & Gam- 
mie 2004; and references therein). In the MHD model 
the magnetorotational instability (MRI; Balbus & Haw- 
ley 1991) drives magnetic turbulence and provides the 
necessary torque to remove angular momentum from the 
gas and drive the inflow. The MHD model is arguably 
the most realistic model for disk accretion of magnetized 
plasma onto a black hole. The standard thin disk model 
provides a simple, analytic, limiting case that is useful 
as a point of comparison. 

In a standard thin disk corotating with the black hole, 
the energy and angular momentum per unit rest mass 
accreted by a black hole are the energy and angular mo- 
mentum of a unit mass at the ISCO, immediately prior 
to its rapid plunge and capture by the hole: 

^2 



-ElSCO : 



2Mr„ 



ay/ Mr u 



3M h r„ 



ISCO : 



2aVM^I) 1/2 
y/Wi~{r 2 ms - 2ay/Wj~ + a 2 ) 

iM h r ms + 2a^Ml~Y/^ 



(9) 
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where the ISCO radius r ms is given by 

r ms = M{3 + Z 2 - [(3 - Zi)(3 + Z l + 2Z2)] 1 / 2 }, (10) 
where 



2 \ 1/3 



1 



M 



and 



1/3 



1/2 



1 



M 



1/3' 
(11) 

(12) 



(see, e.g., Shapiro k Teukolsky 1983, eqns 12.7.17-12.7.18 
and 12.7.24). The mass accretion efficiency and spin evo- 
lution parameters corresponding to the thin disk model 
are then given by 

l--Eisco, (13) 



s — hsco — 2— £isco 



(standard thin disk). (14) 



The MHD disk accretion model of Gammie, Shapiro 
k McKinney (2004) and McKinney & Gammie (2004) is 
based on a fully relativistic, axisymmetric simulation of 
a nonradiative, magnetized plasma onto a Kerr-Schild 
black hole within the MHD approximation. The ini- 
tial gas configuration is a torus with an inner radius at 
r/M = 6; in the absence of a magnetic field the torus is 
constructed to be in equilibrium about the hole (Fish- 
bone & Moncrief 1976). The torus in threaded with 
a poloidal magnetic field initially and evolves with an 
adiabatic equation of state (EOS) with an adiabatic in- 
dex F = 4/3 (to model a radiation-dominated, inner- 
disk EOS). The source of viscosity is MHD turbulence 
driven by the MRI instability. The simulations are per- 
formed for various black hole spin parameters a/M, hold- 
ing the value of the spin parameter fixed during the sim- 
ulation. The evolution proceeds for many dynamical 
timescales M, until a crude steady-state, with fluctua- 
tions, is achieved. 

The results of the numerical simulations suggest that 
in steady-state the radiation efficiency parameter cm as a 
function of a/M is remarkably close to the function char- 
acterizing the standard thin disk fean. IT3*|) . even though 
there is no sharp transition in the surface density at or 
near the ISCO. However, the spin evolution parameter 
s(a/M) is different and can be represented reasonably 
well by the least squares linear fit 



s = 3.14-3.30- 



M 



(MHD disk), 



(15) 



(see Table 2 in McKinney & Gammie 2004). The nu- 
merical simulations demonstrate that the above parame- 
ters describing steady-state, MHD accretion-disk behav- 
ior are not particularly sensitive to the initial conditions 
in the disk (e.g. the initial B— field). As McKinney & 
Gammie (2004) discuss, the key results are also quite 
comparable to those found by De Villiers et al. (2004), 
who used a different numerical method and took the adi- 
abatic index of the gas to be T = 5/3 instead of T = 4/3. 
We therefore model a relativistic MHD accretion disk by 
adopting eqns (|13[) and l|15(l in our evolution equations. 

The s vs. a/M curve for the MHD disk is roughly par- 
allel to, but somewhat below, the curve for the standard 
disk (see Fig 4 in Gammie, Shapiro & McKinney 2004 for 



a comparison). In particular, the parameter s never falls 
to zero for a standard thin disk until the hole is max- 
imally rotating at a/M = 1, while for an MHD disk s 
crosses zero at a/M w 0.95. The crucial physical conse- 
quence is that steady accretion via a standard thin disk 
always causes the central black hole to spin up until it is 
maximally rotating (Bardeen 1970), while accretion via 
an MHD disk drives the black hole to spin equilibrium 
at a/M k, 0.95. While the value of a/M at equilibrium 
is not determined precisely and may depend on details of 
the flow geometry (and the presence of radiative cooling), 
its departure from maximal is robust and is typically 
> 0.05 below unity. The difference between a/M = 1 
and a/M = 0.95 is quite substantial physically: for ex- 
ample, the efficiency paramter is em = 1 — 1/3 1 / 2 = 0.42 
for a/M — 1, while it is considerably smaller, cm = 0.19, 
for a/M — 0.95. As we will see below, this difference 
will have significant consequences for the cosmological 
growth of a black hole from an initial seed. For compari- 
son we note that the radiative efficiency for disk accretion 
onto a nonrotating Schwarzschild black hole (a/M = 0) 
is e M = 1- (8/9) 1 / 2 = 0.057. 

If we assume that the supply of gas remains sufficiently 
copious in the vicinity of the black hole, then the ac- 
cretion luminosity is likely to be Eddington-limited and 
nearly constant, with ej, i=a 1. Accretion models charac- 
terized by super-Eddington luminosities with cl > 1 are 
possible theoretically (Ruszkowski & Begeleman 2003). 
However, the estimated Eddington ratios of the bolo- 
metric luminosities of the broad-line quasars in a Sloan 
Digital Sky Survey sample of 12,698 quasars in the red- 
shift interval 0.1 < z < 2.1 support the value cl ~ 1 as a 
physical upper limit (McLure & Dunlop 2004). (Here the 
mass of the central black hole is estimated from a virial 
relation, assuming that the gas motions in the broad-line 
region of the quasars, measured by H/3 and Mgll emis- 
sion lines, are virialized.) The maximum allowed value of 
cl establishes the maximum growth rate of a black hole 
from accretion (see eqn. |SJ) and, as a result, the most 
robust (least stringent) constraints on any gas accretion 
scenario for the growth of a supermassive black hole from 
a smaller intial seed. To establish these constraints we 
therefore shall set — 1 in many of our numerical esti- 
mates below. 

As a point of reference, and for later application, it is 
useful to integrate eqn. JSJ) assuming that the mass and 
luminosity efficiencies both remain constant with time, 
yielding 

M(t)/M&)=exp '^(l-e M )(t- ti y 



{^m^l constant), 
(16) 



where ti is the initial time at which the black hole has 
a mass Mi. Note that the right hand side of eqn. l|16fl 
is independent of black hole mass. This fact makes it 
possible, under certain plausible conditions that we shall 
specify, to disentangle and track separately the amplifi- 
cation of black hole mass by accretion from the ampli- 
fication by discrete mergers. Let M n (t) be the mass of 
the black hole at time t following its n th merger with 
another hole at time t n , where t n < t < t n+ i. Assume 
that the duration of a merger, as well as the time re- 
quired for accretion to drive the merged remnant to spin 
equilibrium (see ean. 1231 below), are both much shorter 
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than the time interval between mergers, and that the 
hole continues to accrete steadily throughout this inter- 
val. Note that black hole mergers can completely eject 
black holes from halo centers owing to gravitational wave 
recoil, and thereby turn off accretion altogether (see dis- 
cussions of black hole recoil in halos in, e.g., Hut & Rees 
1992; Merritt et al. 2004; Madau & Quataert 2004). 
However, incorporating the most recent recoil calcula- 
tions (Favata, Hughes & Holz 2004) into simple models 
of dark halo mergers, Yoo & Miralda-Escude (2004) con- 
clude that the kick velocities are not sufficiently large 
to impede black hole growth significantly (cf. Haiman 
2004). Note also that a major merger between two black 
holes of comparable mass may change both the magni- 
tude and direction of the spin of the resulting black hole 
remnant. Following such a merger, the orientation of 
the black hole spin may not be aligned with the orienta- 
tion of the asymptotic gaseous disk at radii r ^> 100M 
outside the hole. (Rees 1978 and Natarajan & Pringle 
1998 point out that the black hole exerts a torque on the 
asymptotic disk, which also exerts a torque back on the 
hole, eventually forcing their mutual alignment, but on 
a timescale that is still uncertain). Moreover, the ori- 
entation of the asymptotic disk will likely fluctuate in 
time, due to the redistribution of gas following mergers 
of dark halo cores, galaxy mergers, and the tidal disrup- 
tions of passing stars by the central hole. However, near 
the black hole, at radii from r ~ M to r ~ 20M, where 
the bulk of the disk's gravitational energy is released and 
the hole-disk interactions are strong, the hole's gravito- 
magnetic field will exert a force on the disk that, when 
combined with viscous forces and magnetic fields, will 
drive the disk down into the hole's equatorial plane, (the 
"Bardeen-Petterson effect"; Bardeen & Petterson 1975). 
This effect will maintain the alignment between the axis 
of the inner disk and the spin axis of the hole. Accretion 
will subsequently drive the hole to spin equilibrium and 
restore eu to its equilibrium value. 

For the situation described above we can take €m to 
be constant and equal to its value at spin equilibrium 
throughout most of the lifetime of the hole. Assume fur- 
ther that is also constant, which should be the case 
if the available gas is sufficiently copious that the accre- 
tion is always Eddington- limited, whereby ~ 1. Let 
f n be the mass amplification of the hole following its n th 
merger with another hole: /„ = M n (t n ) /M n -\(t n ) > 1. 
Then we may use eqn. Ijl6(l to calculate the total mass 
amplification from t{ to tf according to 



M f = M N (t f ) = Mpfti) Mi(ti) Mi(t 2 ) 
Mi M (U) M (U) Mo(h) Mi(ti)"' 
M N -i(t N ) M N {t N ) M N (t f ) 



(17) 



Mjv_i(i w _i) M N -i{t N ) M N {t N ) 
= exp [C(ti - U)] h exp [C(t 2 - h)} ■ ■ • (18) 

• • - exp [C(t N - tjv_i)] f N exp [C(t f - t N )] 
= hh ■■■ In exp [C(t f - U)] , 

where C — cl(1 — eu) / (cmt) and is essentially constant. 
Comparing eqns. (|16[) and (|17fl reveals that the net mass 
amplification due to accretion can be treated as a single 
multiplicative factor that is independent of the net am- 
plification factor due to mergers, / = /1/2 • • • In- Thus, 
for the simple scenario envisioned here, eqn. (|16(l can 



be applied to determine the growth of a black hole by 
accretion, even when steady growth by accretion is in- 
terrupted and augmented by discrete, stochastic black 
hole mergers. 

2.2. Cosmological Model 

To relate the time parameter t appearing in the evo- 
lution equations to an observable time-like quantity like 
the redshift z, we adopt the concordance ACDM, spa- 
tially flat (k = 0), cosmological model. All the free pa- 
rameters in this model that we shall need for our compu- 
tations have been measured by now, so that the model 
is uniquely specified. 

The basic evolution equation for the Fricdmann- 
Robertson- Walker expansion parameter a(t) is given by 
2 



H 2 = 



- =m 



n 



fooY 
V a J 



(19) 



where H (t) is Hubble's constant and the normalized mass 
density parameter fij^ and cosmological constant param- 
eter $7° satisfy the relation 

C + n° A = i. (20) 

In the above equations the sub(super)script "0" denotes 
the value of a quantity at the current epoch, z = Q. 
Recalling that a/a = 1 + z and substituting eqn. I|20|l 
into eqn. I|19|) . we can integrate eqn. I|19|) to obtain t — 
t(z): 



t(z) 



3ff (l-^„) 1 / 2 



sinh 



1 



(1 + zf 



1/2 



(21) 

To evaluate t(z) appearing above we adopt the WMAP 
values for the cosmological parameters: fi^n ~ 0.27 and 
H = lOQ/i km/s/Mpc with h w 0.71 (Bennett et al. 
2003; Spergel et al. 2003). 

To evaluate the molecular weight /i e appearing in 
eqns. © and © we assume that the accreting plasma is 
a zero-metallicity, primoridial gas, for which 

* = Y~Y/2' (22) 

We take the primordial helium abundance Y to be Y f=a 
0.25 (Cyburt, Fields & Olive 2003). 

For the concordance model, the age of the universe is 
T = t(0) = 13.7 Gyr, while the age at redshift z = 6.43, 
the highest known quasar redshift, (SDSS 1148+525; Fan 
et al. 2003), is only £qso = 0.87 Gyr. Hence £qso is the 
upper limit to the time available for accretion to occur 
onto the initial seed black hole that powers this quasar. 
In fact, assuming that the black hole seed forms from 
the collapse of a first generation, Pop III star at redshift 
z < 40 and t > f(40) = 0.067 Gyr, the available time 
for accretion is reduced to f accrete ^0.80 Gyr. We note 
that the stellar evolution (hydrogen burning) lifetime of 
a massive Pop III progenitor, t cvo ~ 0.003 Gyr (Wag- 
oner 1969; Onukai & Palla 2003), is much smaller than 
t(40), so that the delay between stellar formation and 
collapse is of little consequence for determining the to- 
tal time available for accretion growth. Most important, 
the exponential accretion growth timescale is given by 
eqns. © and JUJ to be 
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which is considerably smaller than t accrete- It is the avail- 
ability of many exponential growth timescales from the 
time of black hole seed formation to the birth of a quasar 
that makes it possible for the black hole to grow from 
stellar to supermassive size by gas accretion (aided by 
mergers) in the early universe. 

3. BLACK HOLE GROWTH AND SPINUP 

Here we integrate the coupled mass and spin evolution 
equations JSJ and © to study black hole growth and 
spin-up by gas accretion. In Fig. 1 we show the increase 
in mass by accretion at the Eddington limit (e^ = 1) as 
a function of time. In (a) the accretion disk is a stan- 
dard thin disk with a "no torque boundary condition" 
at the ISCO (eqn. O; in (b) the disk is an MHD disk 
as calculated by McKinney & Gammie (2004) (eqn. 1 15(1 . 
We consider two different initial values for the black hole 
spin parameter, a/M. The solid lines show the case in 
which the initial black hole is nonrotating with a/M = 0. 
The dotted lines show the case in the which the black 
hole to be spinning with a/M = 0.75. The latter is the 
value calculated for a black hole formed from the catas- 
trophic collapse of a massive, radiation-dominated star 
spinning uniformly at the mass-shedding limit that has 
evolved to the onset of relativistic radial instability just 
prior to collapse (Shibata & Shapiro 2002; Shapiro & 
Shibata 2002; Shapiro 2004). For each disk model, the 
solid dots show the black hole growth that would oc- 
cur if the accretion were maintained at the asymptotic 
black hole spin and efficiency from the beginning. Time 
is plotted in units of r given by eqns. I© and (|22|l for 
a zero-metallicity, cosmological abundance of H and He. 
The total duration plotted corresponds to the age dif- 
ference between redshift z% = 40 (the earliest plausible 
birthdate of a black hole from the collapse of a Pop III 
star) and redshift Zf = 6.43 (the highest known quasar 
redshift, corresponding to SDSS 1148+5251; Fan et al. 
2003) in the adopted ACDM cosmology. 

The key difference between the two disk models is that 
the standard thin disk drives the black hole to max- 
imal spin and efficiency (a/M = 1 and cm = 0.42; 
Bardeen 1970) while the MHD disk drives the black hole 
to spin equilibrium at lower values (a/M « 0.95 and 
ejvf ps 0.19). A moderately lower radiation efficiency em 
for Eddington-limited accretion results in substantially 
larger black hole growth. 

As is evident from the figure, apart from an initial 
transient lasting ~ O.lr ~ r grw th (see figure inset), the 
asymptotic evolution and final black hole mass amplifica- 
tion factor is little affected by the initial black hole spin. 

In Fig. 2 we show the spin evolution during Eddington- 
limited accretion for the cases shown in Fig. 1. The large 
figure plots the spin parameter versus time, while the in- 
set plots the spin parameter as a function of black hole 
mass amplification. After the initial transient lasting ~ 
O.lr, during which time the black hole grows by a factor 
of ~ 2, the black hole spin and efficiency approach their 
asymptotic values. An analytic integration of eqn. JHJ 
with eqn. Ijl5|l , crudely holding cm fixed, yields an ex- 
ponential spin-up timescale r sp i n sa (1 — £M)T grw th/3.30, 
which explains the rapid spin-up rate. This result also 
explains why the asymptotic black hole evolution and 
mass amplification are little affected by the initial black 
hole spin and depend only on the equilibrium spin rate 




Fig. 1. — Black hole growth by accretion at the Eddington limit. 
In (a) the accretion disk is a standard thin disk with a "no torque 
boundary condition" at the ISCO, ignoring photon recapture; in 
(b) the disk is an MHD disk as calculated by McKinney & Gam- 
mie (2004). The solid lines show evolution starting from a/M = 0, 
while the dotted lines show evolution starting from a/M = 0.75. 
For each disk model, the solid dots show the mass amplification 
that would occur if the accretion were maintained at the asymp- 
totic, equilibrium spin and radiation efficiency from the beginning. 
Time is plotted in units of r = Mc 2 /Le = 0A5/fi e Gyr. The 
total duration shown corresponds to the age difference between 
redshift zi = 40 (the earliest plausible birthdate of a black hole 
from the collapse of a Pop III star) and redshift Zf = 6.43 (the 
highest known quasar redshift, corresponding to 1148+5251) in a 
concordance ACDM cosmology. 



of the accreting Kerr black hole. Thus, for all but the 
initial transient, it is adequate to adopt the equilibrium 
spin rate and use eqn. I|16fl with the corresponding con- 
stant efficiency eA/[(a/M) cq ] to calculate the cosmolog- 
ical black hole mass amplification in lieu of a detailed 
integration of coupled evolution equations. 

4. COSMOLOGICAL IMPLICATIONS 

Next we consider some of the cosmological implications 
of the results of the previous section for the growth of 
SMBHs in the early universe. 

In Fig. 3 we employ eqns. H16[) and (|21|l to study 
the black hole mass amplification factor Mf/Mi as a 
function of the redshift z% of the initial seed black hole. 
Here we plot the final amplification factor achieved by 
redshift Zf = 6.43, the highest known quasar redshift. 
Each solid curve is labelled by the adopted constant ra- 
diation efficiency, e^] the luminosity is assumed to be 
the Eddington value (el = 1). The horizontal dashed 
lines indicate the range of amplification factors required 
for accretion alone to grow a supermassive black hole 
of mass 10 9 M Q from an initial seed black hole of mass 
100 < M/Mq < 600 formed from the collapse of a Pop 
III star. A mass of 10 9 M Q is the value inferred for typi- 
cal quasars, including 1148+5251 (Fan et al. 2003). The 
lower the mass of the initial seed, the larger is the re- 
quired mass amplification. The horizontal dotted lines 
indicate the required range of accretion-driven mass am- 
plification assuming that mergers assist the growth pro- 
cess by accounting for an amplification of / ~ 10 4 in 
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Fig. 2. — Spin evolution during Eddington-limitcd accretion for 
the cases shown in Fig 1. Curves are labelled as in that figure. The 
large figure plots the spin parameter versus time, while the inset 
plots the spin parameter versus the black hole mass amplification 
factor. 
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Fig. 3. — Black hole accretion mass amplification Mr /Mi ver- 
sus redshift Zi of the initial seed. Here we plot the amplification 
achieved by redshift Zf = 6.43, the highest known quasar red- 
shift, corresponding to 1148+5251. Each solid curve is labelled 
by the adopted constant radiation efficiency, ejvf; the luminosity 
is assumed to be the Eddington value (e^ = 1). The horizontal 
dashed lines bracket the range of amplification required for accre- 
tion alone to grow a seed black hole of mass 100 < M/Mq < 600 
formed from the collapse of a Pop III star to 10 9 Mq. The horizon- 
tal dotted lines bracket the required accretion amplification range 
assuming that mergers account for a growth of 10 4 in black hole 
mass, the remainder being by gas accretion. 



black hole mass, with the remaining amplication pro- 
vided by gas accretion (Yoo & Miralda-Escude 2004). 

Some important cosmological implications can be in- 
ferred from Fig 3. Consider a seed black hole that forms 
sometime after redshift Zj < 40, by which time the ear- 
liest stars have formed and collapsed. In the absence of 
mergers, steady accretion cannot by itself achieve the 



required growth to explain quasars at Zf = 6.43 un- 
less the efficiency satisfies eu < 0.13. For disk accre- 
tion, this constraint requires the steady-state black hole 
spin parameter to be below a/M w 0.83. This value is 
below (1) the maximal black hole spin a/M = 1, the 
asymptotic equilibrium value for a standard thin disk; 
(2) a/M = 0.998 (for which e M ~ 0.32), the equilibrium 
value of a standard thin disk accounting for the recap- 
ture of some of the emitted photons by the black hole 
(Thorne 1974); and (3) a/M w 0.95, the equilibrium 
value of a typical MHD disk. Therefore, it is likely that 
mergers are required to assist accretion to achieve black 
hole growth to supermassive size by zt = 6.43. As may 
be inferred from the figure, this conclusion holds even in 
the (unlikely) event that the black hole seed forms much 
earlier than Zi w 40. 

Monte Carlo simulations by Yoo & Miralda-Escude 
(2004) of hierarchical CDM halo mergers, accompanied 
by mergers of their central black holes, suggest that black 
hole mass amplification factors of ~ 10 4 are achieved 
via mergers by Zf = 6.43. This result implies that 
/ = /1/2 • • • In ~ 10 4 in eqn. JTZ). While other sim- 
ulations predict lower growth factors due to mergers, we 
adopt 10 4 to be conservative, noting that our conclu- 
sions below are strengthened if the smaller factors turn 
out to be correct. As eqn. i|17|) shows, eqn. (|16|) can be 
used to estimate the cumulative mass amplication by gas 
accretion, even when steady growth by accretion is inter- 
rupted and enhanced by stochastic black hole mergers. 
Minor mergers of massive black holes with smaller holes 
tend to spin-down the massive hole (Hughes & Bland- 
ford 2003; Gammie, Shapiro & McKinney 2004). Ma- 
jor mergers of two black holes of comparable mass fol- 
lowing binary inspiral drive the the merged remnant to 
a/m > 0.7. [See, e.g., Baumgarte & Shapiro 2003 for a 
review of binary black hole coalescence, and references; 
see Gammie, Shapiro & McKinney 2004 for a general dis- 
cussion of black hole spin-up and spin-down mechanisms, 
and references]. In either case, after a short transient 
epoch, accretion will drive the merged remnant to the 
disk accretion equilibrium spin rate and corresponding 
mass efficiency (see Section 2), as assumed in employing 
eqn. (JTBJ). 

Fig 3. indicates that for black hole seeds that arise 
from collapsed Pop III stars, merger-assisted mass am- 
plification to SMBH status is easily achieved for typical 
MHD accretion disks, marginally possible for a standard 
thin disk in spin equilibrium accounting for photon re- 
capture, but not possible for a standard thin disk that 
drives the black hole to maximal rotation. However, if 
the black hole seed is less than 600Mq ; the standard thin 
disk appears to be ruled out: the required amplification 
places the lower dotted line in the plot above the curve 
for ejvf = 0.32, corresponding to a/M = 0.998, the value 
for a standard thin disk in spin equilibrium, accounting 
for photon recapture. These conclusions again also hold 
even if the black seed forms much earlier than zi k, 40, 
while they may be tightened if the seed forms later. In 
fact, it may be likely that the seed forms later, at Zi < 40, 
given that even 4 — a peaks in the density perturbation 
spectrum for the progenitor halo of SDSS 1148+5251 do 
not collapse until z ~ 30 in the ACDM concordance cos- 
mology (see, e.g. Figure 5 in Barkana & Loeb 2001). 
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Moreover, the potential wells of the earliest halos are 
quite shallow (~ 1 km/s) and may not be able to retain 
enough gas to form stars. Nevertheless, the effect of al- 
tering the date of birth of the BH seed is not very great 
unless Zi < 20 — 25, as is evident from Fig. 3. 

It is significant that the range of equilibrium accretion 
disk radiation efficiences required to achieve the neces- 
sary growth of a black hole seed to supermassive size by 
Zf = 6.43 is consistent with the values inferred observa- 
tionally for R, the ratio of the QSO plus AGN luminosity 
density to the mass density of SMBHs in nearby galaxies: 
e M > R ~ 0.1 - 0.2 (Soltan 1982; Yu & Tremaine 2002; 
Elvis et al. 2002). This consistency supports the notion 
that accretion plays an important role in SMBH growth 
and is responsible for the acquisition of the bulk of the 
final mass of a SMBH. The range of cm inferred from 
R favors accretion disk models that drive the black hole 
to spin equilibrium in the range 0.7 < a/M < 0.95, well 
below maximal rotation and consistent with the values 
found in recent simulations of relativistic MHD accretion 
disks. 

How are our conclusions altered if (when) a quasar is 
discovered at a higher redshift, zqso > 6.43? We antic- 
ipate this possibility in Fig. 4, where we solve eqn. l|16fl 
to plot the mass efficiency ejf required to build a super- 
massive black hole from a seed as a function of the host 
quasar redshift, zqso- Here we again assume that a seed 
black hole of mass 100 < M/M Q < 600, formed from 
the collapse of a Pop III star at Zi 3> zqso, grows by 
zqso to 10 9 M Q . While accreting, the black hole again 
is assumed to radiate at the Eddington luminosity, with 
6l = 1. We set ti — 0, Zi — oo, so that the plotted radia- 
tion efficiency represents an upper limit to em ; increasing 
ti and lowering Zi decreases the required efficiency (but 
the decrease is small if, for example, we set Zi > 40). The 
dashed lines indicate the range of efficiencies required for 
gas accretion alone to achieve the necessary growth. The 
dotted lines indicate the range required for gas accretion 
assuming that mergers assist the growth and account for 
an amplification of / = 10 4 in black hole mass. The hori- 
zontal solid lines bracket the range of efficiencies inferred 
observationally for R. 

In the absence of mergers, the upper limit to the ra- 
diation efficiency required to build a supermassive black 
hole by zqso ~ 6.43 is em ~ 0.14, consistent with obser- 
vational constraints implied by R, but inconsistent with 
the much higher values expected for a standard thin disk 
driven to spin equilibrium, and marginally inconsistent 
for an MHD accretion disk. The inconsistencies become 
worse as zqso > 6.43. In fact, the upper limit to em 
falls below 0.1 for zqso ^ 8, which would force cm to lie 
below the observationally inferred range for R and would 
therefore be difficult to understand. 

When mergers are included, the upper limit to the ef- 
ficiency required to build a SMBH by zqso ~ 6.43 in- 
creases to €m ~ 0.30 (i.e., any black hole seed born at a 
finite redshift, oo > z\ > zqsOi must accrete at a lower 
efficiency than 0.30 to reach 10 9 M Q ). This upper limit is 
roughly consistent with the observational constraint on 
cm > R and with theoretical values for accretion from 
an MHD disk in spin equilibrium, but only marginally 
consistent for accretion from a standard thin disk in spin 
equilibrium accounting for photon recapture, and incon- 
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Fig. 4. — The radiation efficiency ejvf required to build a SMBH 
versus the redshift zqso of the host quasar. The seed black hole 
forms with a mass 100 < M/Mq < 600 from the collapse of a 
Pop III star at Zi S> zqso an d grows by zqso t° 10 9 Mq. Accre- 
tion occurs at the Eddington luminosity. The dashed lines bracket 
the range of efficiencies required for accretion alone to achieve the 
necessary growth. The dotted lines bracket the range required for 
accretion assuming that mergers account for a growth of 10 4 in 
black hole mass. The horizontal solid lines bracket the range of 
efficiencies inferred from R, the observed ratio of the QSO plus 
AGN luminosity density to the mass density of local SMBHs. 



sistent for accretion from a standard thin disk that drives 
the black hole to maximal rotation. Should a quasar be 
discovered at zqso > 6.43, it would appear that accre- 
tion from a standard thin disk will be ruled out: the upper 
limit will fall below em = 0.32, the value for a standard 
thin disk in spin equilibrium accounting for photon re- 
capture. Should a quasar be discovered at zqso ^ 10, em 
would fall below 0.19 and the results would be difficult to 
reconcile with accretion from a typical MHD disk as mod- 
eled in recent simulations. These critical values of zqso 
are all smaller if the initial black hole seed is smaller than 
600Mq, as the top curve, which sets the limit, is lowered 
in the figure. Finally, should a quasar be discoverd (per- 
haps unexpectedly) at zqso ^ 18, the upper limit to em 
would drop below the observationally inferred value 0.1 
for all Zi > zqso- This drop would not be easily ex- 
plained by any of these models of accretion, and might 
require super-Eddington accretion with el > 1 to keep 
em > 0.1. Alternatively, the value of R, observed for 
local host galaxies, might be significantly smaller at high 
redshift, which would relax the constraint on em- But 
achieving the inferred lower radiation efficiencies might 
then require ADAF disks or spherical accretion. 

5. SUMMARY 

We have explored the evolution of black hole mass and 
spin by gas accretion in the early universe. We have illus- 
trated how for Eddington-limited accretion, the growth 
of a SMBH depends sensitively on the radiative efficiency, 
em ■ For disk accretion, the mean efficiency is determined 
by the equilibrium black hole spin, which in turn depends 
on the the torques acting on the gas near the black hole 
horizon. We have explored the consequences of the as- 
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sumptions that seed black holes are the remnants of col- 
lapsed Pop III stars that may form as early as z < 40 
and can grow to ~ 10 9 M e by zqso = 6.43, the highest 
redshift discoverd to date, corresponding to QSO SDS 
1148+5251. Allowing for growth both by accretion and 
mergers, simple theory suggests that the required mass 
amplification is possible provided the radiation efficiency 
satisfies eu < 0.2, with the upper limit decreasing should 
a quasar be discoverd at higher redshift zqso > 6.43. 
The inferred efficiency is consistent with the observed 
ratio R of the QSO plus AGN radiation density to the 
mass density of SMBHs in local galaxies, which suggests 
that an appreciable fraction of the mass of the final black 
hole is acquired by disk accretion, rather than by merg- 
ers. The inferred efficiency favors MHD accretion disk 
models that exert non-zero torques on the gas at the in- 
ner edge of the disk. These disks ultimately drive the 
black hole to spin equilibrium at a/M ~ 0.95, substan- 
tially below the maximum spin allowed for a Kerr black 
hole. 

We have made a number of simplying assumptions in 
an effort to provide a concrete calculation that illustrates 
how the SMBH initial seed, its cosmological growth, and 
accretion disk models may all be constrained by the 
existence of quasars at high redshift. Our conclusions 
are tentative, as there exist many caveats. For exam- 
ple, if the progenitors of black hole seeds are SMSs with 
Al ^> 10 3 M Q instead of Pop III stars, many of the con- 
straints on the radiation efficiency cm, and the associ- 
ated range of expected black hole spins, may have to 
be relaxed. On the other hand, if the progenitors have 
masses < IOOMq the constraints favoring low efficien- 
cies £m 0.2 characterizing MHD accretion disks in 
spin equilibrium are strengthened. If the accretion lumi- 
nosity is super-Eddington with ez, > 1, then the upper 
limit for cm increases (see eqn. 1161) . relaxing the con- 
straints which favor an MHD disk over a standard thin 
disk. Super-Eddington accretion is possible theoretically 
(Ruszkowski & Begelman 2003). However, the a Sloan 
Digital Sky Survey of 12,698 broad-line quasars in the 
redshift interval 0.1 < z < 2.1 supports the Eddington 
value as a physical limit (McLure & Dunlop 2004). More- 
over, given the observed value for R, a higher limit for 
cm would imply that a large fraction of the emitted ra- 
diation lies outside the optical bandwidths included in 
the determination of i?, which seems unlikely. Finally, 
mergers may lead to black hole mass amplification fac- 
tors smaller than / ~ 10 , the typical value found in 
the Monte-Carlo calculations of Yoo & Miralde-Escude, 
and the value assumed here. A smaller merger amplifica- 
tion factor requires larger accretion-driven amplification 
in order for a black hole seed to reach SMBH size by 
z qso — 6.43 (see eqn. IT?}) . This would strengthen the 
constraints set by SMBH accretion growth on the upper 
limit to ej\/, although the efficiency must still be consis- 
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tent with R if the bolometric accretion luminosity resides 
mostly in the optical wavebands included in R and the 
value of R measured for local galaxies also applies at 
higher redshift. Of course, the peak contribution to the 
value of R arises from the average behavior of accreting 
BHs with z ~ 3; this ratio does not necessarily apply 
to individual objects like SDSS 1148+5251, nor need it 
apply to SMBHs at z 6. 

We have assumed that the mass of SDSS 1148+5251 
is ~ 10 9 Af Q ; a lower value would relax many of our con- 
straints, which a higher value would strengthen them. If 
the flux from this source were amplified by gravitational 
lensing, or beaming, then a lower mass estimate would be 
appropriate. However, no multiple images have been seen 
(Richards et al. 2004) and it has been shown that high 
amplification without at least two images is very improb- 
able (Keeton et al. 2004). Strong beaming also seems 
unlikely since it would reduce the line/continuum ratio 
(Haiman & Cen 2002), which is not observed (Willott et 
al. 2003). In fact, assuming that the quasar emits at 
the Eddington luminosity gives a mass of 4.6 x 10 9 M© 
(Fan et al. 2003; Haiman 2004), and this higher value 
strengthens our conclusions somewhat. 

Isolating the accretion growth of a seed black hole in 
the early universe from the hole's full dynamical history 
and environment does not allow us to account for other 
important correlations that provide clues to the forma- 
tion of SMBHs. Such correlations include the SMBH 
mass vs. bulge luminosity relation, Mbh oc ibuigo (Kor- 
mendy & Richstone 1995), and the SMBH mass vs. ve- 
locity dispersion relation, Mbh °c (Gebhardt et al. 
2000; Ferrarese & Merritt 2000; Tremaine et al. 2002), 
inferred for nearby host galaxies. Only by performing de- 
tailed simulations that track the formation and growth 
of SMBHs in a cosmological setting governed by hier- 
archical halo mergers, black hole mergers, gas settling, 
star formation, and feedback can these correlations be 
reliably reproduced. We look forward to the next gener- 
ation of simulations that incorporate the recent results 
of relativistic MHD accretion onto black holes, since, as 
demonstrated here, the outcome of these global simula- 
tions may depend sensitively on the local physics of such 
accretion flows. 
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